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A two-dimensional version of the shear-transformation-zone (STZ) theory by Falk and Langer is 
explored numerically. Two different geometries are used to simulate uniaxial tension experiments 
where materials are subjected to constant strain rates. In the first setup, a rectangular specimen is 
given an imperceptible indentation, allowing it to neck at the center. The dynamics are explored 
systematically by varying both the straining capability of the STZs (eo) and the external strain 
rate. Higher values of eo increase the plastic flow and result in sharper necks. Decreased values of 
the external strain rate result in the formation of narrow shear bands, consistent with the absence 
of thermal relaxation mechanisms in the model. In the second configuration, the equivalent of pre- 
annealed materials are explored. Here, the sample is initially square and the edges are made rough 
in order to encourage the formation of shear bands. With respect to eo and the external strain rate, 
the results show trends similar to those in the necking simulations. The pre-annealing was modeled 
by using a low initial density of STZs. This had most effect when eo was large, contributing to the 
localization of the strain and making the material appear more brittle. 
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I. INTRODUCTION 

The goal of the continuum shear-transformation-zone 
(STZ) theory is to describe plastic deformation in amor- 
phous solids on a mesoscopic scale, averaging out some of 
the microscopic, discrete details. The plastic deformation 
is described in terms of flow rates, similar in approach to 
that of the Navier-Stokes model pj although the scalar 
pressure has been replaced with a stress tensor. 

The STZ theory was constructed by Falk and Langer 
H H II H H IS S H and originated with the assump- 
tion that plastic deformation is limited to, and defined 
as, the non-affine transformation of particles in localized 
areas or zones. This was based on similar ideas made 
by Argon, Spaepen, and others who described creep in 
metallic glasses in terms of local, molecular transitions 
or rearrangements [3 [3 13 13 13 El This, in turn, 
grew out of theories by Turnbull, Cohen, and others, who 
suggested that the observed behavior in the amorphous 
metals could be described by lin king the transition rates 
to local free- volume fluctuations 0, 0, 0, . 

This paper, continuing the work reported in |19| . 
presents a two-dimensional model that describes elastic 
and plastic dynamics in an amorphous solid; the STZ 
theory supplies the model with a plastic-flow description. 
The main contribution of this paper is the exploration of 
the STZ theory in a spatially extended system, focus- 
ing on features introduced by having an inhomogeneous 
geometry. Section [H] presents a tensorial version of the 
STZ theory that was developed in 20] and based on both 
the scalar model developed earlier by Falk and Langer, 



as well as some of Falk's initial ideas of how to write the 
theory in a two-dimensional setting 0. 

The remaining part, which forms the core of the paper, 
reviews two sets of numerical uniaxial-tension constant- 
strain-rate simulations of the two-dimensional, tensorial 
STZ model. First, Section ILTT1 gives a brief discussion of 
the implementation. Then the first set is presented in 
Section llVl exploring how the plastic flow affected the 
dynamics during necking. The second set is described in 
Section[V] relating the pre-annealing of amorphous mate- 
rials to increased strain localization and brittle behavior. 
The paper is concluded in Section I VII 



II. THE TENSORIAL STZ MODEL 

The numerical simulations described in this paper uses 
the so-called quasi-linear tensorial STZ model in two di- 
mensions, given below. When combined with boundary 
conditions, it can be used to describe both the elastic and 
plastic dynamics of an amorphous solid in a spatially ex- 
tended geometry. The equations are 
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The variables here are the pressure p, the deviatoric 
stress sy, the velocity V{, the relative density of STZs 
A, and the STZ alignment A ir D\f and Df- are total 
and plastic rates of deformation, respectively. The tilde 
on the former indicates that only the deviatoric part of 
the tensor is being used (that is, the trace has been sub- 
tracted from it), r is proportional to both the rate at 
which STZs arc created and annihilated, as well as to 
the rate at which energy is dissipated through plastic 
deformations. The invariant of the rate of plastic defor- 
mation is defined as D pl = J The parameters 
in the equations, which remain constant throughout a 
simulation, are the shear modulus /i, the two-dimensional 
Poisson ratio V2, the mass density p, and the straining ca- 
pability of the STZs eo- The temporal derivatives D/Dt 
and T>/T>t both include an advective term, and the latter 
also takes rotation into account. The complete deriva- 
tion of these exact equations can be found in J2jj; see 

also uimnm. 

Concerning expressions for the stored energy, it was 
found in an earlier publication |2l| that the plastic dissi- 
pation is given by [46( 



Q = 2e Ar , 
and the plastic energy density 
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In addition, the elastic and kinetic energy densities are 
given by 
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One possible interpretation of how plastic deformation 
might be described by the STZ model, involves the flip- 
ping of STZs. A deformation caused by some STZs flip- 
ping, will encourage the internal stress distribution to 
change, which in turn will flip even more STZs as well 
as create new ones and annihilate existing ones. It takes 
energy to flip an STZ. The work done on it could even 
be considered reversible, if one could make the STZ flip 



back again. Then again, the STZ could make it back 
to its original position without releasing any energy by 
being annihilated and recreated, corresponding to perma- 
nent deformation. Since, theoretically, one could regain 
the energy in the flipped STZ if one could make all the 
STZ flip back again without any of them being annihi- 
lated, a reasonable interpretation of stored plastic energy 
ijj pl could be the energy stored in the flipped STZs. In 
practice, though, plastic deformation is an irreversible 
process; some of the stored plastic energy would be lost 
if the STZs were to flip back, since this flipping would 
cause further deformations and thus more creation and 
annihilation of STZs. The plastic dissipation Q can be 
interpreted as the energy lost when an STZ is annihilated 
and recreated. 

Simplified, zero-dimensional models using a non- 
tensorial version of the STZ equations have been ex- 
plored in earlier publications U 0, 01 ■ The cur- 
rent manuscript concentrates on the effects of a two- 
dimensional geometry. 



III. NUMERICAL SETUP 

In order to explore the tensorial version of the STZ 
theory given in Eqs. |JT)|. a couple of two-dimensional 
geometrical configurations were implemented and sim- 
ulated numerically. For this purpose, a C-l — I— program 
using finite-difference algorithms on a regular grid with 
a second-order explicit time-stepping scheme to integrate 
the equations was written from scratch. This section will 
give a brief overview over the numerical implementation 
of the model, including the mapping of the variables onto 
a unit square, numerical algorithms and approximations, 
and boundary conditions. 

The two-dimensional STZ theory was implemented nu- 
merically in order to simulate uniaxial tension experi- 
ments where material specimens were strained at con- 
stant rates. The velocity was controlled along the top 
and bottom edges (the "grips" ) , while the left and right 
boundaries were assumed to have no normal stresses and 
were allowed to deform; see Fig. ^ 

Two sets of simulations were performed; in the first set 
the material was made to neck fSection llVf) . while in the 
second set the samples were given rough boundaries to 
encourage the formation of shear bands (Section 0). In 
the necking simulations, a small indentation was applied 
at the center of the material to break the symmetry and 
to encourage the neck to form in the middle. In addi- 
tion, the material was assumed to be symmetric across 
both the x- and j/-axis, so only a quarter of the system 
was modeled numerically. In the simulations with rough 
boundaries, the left and right edges were made stress- 
free and independent of each other, and no symmetry 
assumptions were used. The bottom of the material was 
held fixed while the top was subjected to a constant strain 
rate. 

One aim when introducing the two-dimensional model 



3 



constant 
strainrate 




yA 



X(y,t) X°(y,t) X'(y,t) 



free boundaries (<?„„= 0) 



FIG. 1: The geometrical setup of the two-dimensional simu- 
lations. A constant strain rate was applied to the rigid ends 
at Y(t) by adjusting the velocity. The free boundaries at 
X(y,t), X°(y,t), and X 1 (y,t) had no normal stresses, that 
is, o nT i = 0. In the necking simulations on the left, the ma- 
terial was assumed to be symmetric across the x and y axes; 
thus only the upper right-hand corner of the material was ac- 
tually simulated. In the simulations on the right, the lower 
edge was held fixed. 



was to find out how it differed from the zero-dimensional 
description, in particular with respect to geometrical in- 
homogeneities. After deciding to investigate constant 
strain rate simulations, it seemed sufficient to only al- 
low the sides to deform, while keeping the grips straight. 
It was assumed that the free boundaries could be de- 
scribed by functions that were single valued: X(y,t) for 
the necking simulations, or X°(y,t) and X (y,t) in the 
case of the model without symmetry assumptions. The 
position of the grips was described by a function only 
dependent on time, Y(t). With these assumptions, the 
material was easily mapped onto a unit square with co- 
ordinates Cx S [0,1] and £ y G [0,1]. For the necking 
simulations with the axis symmetries, the mapping was 
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while the simulations with the free boundaries needed the 
slightly more general transformation 
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See Fig. [21 for an illustration. 

The advantage of mapping the coordinates using 
Eqs. I© and Q) was that the calculations could be done 
on a regular grid where the values of the fields were dis- 
cretizcd into equally spaced nodes. The disadvantage was 
that the operators, like the first-order derivatives, were 
more complicated to calculate. 
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FIG. 2: For the necking simulations, the deformed area 
x G [0, X(y,t)], y G [0, Y(t)] was mapped onto the unit 
square £ x G [0,1], G [0,1]. Similarly, the non- 

symmetric simulations (not shown here) mapped the area 
x G [X°(y,t),X 1 (y,t)] } y G [0, Y(t)) onto the same unit 
square. 



All the simulations were done at constant strain rate. 
This was enforced by controlling the y-component of the 
velocity, v y , at the grips (the top and bottom boundaries 
in Fig. For example, in the necking simulations the 
grip velocity was given by 



l [x,y = ±Y(t),t] = ±Y(t)D* 



(8) 



where D tot was the constant strain rate. Since the veloc- 
ity was applied to all the nodes along the top and bot- 
tom, these edges remained straight. The position of the 
edge as a function of time, using that v y [x, y — Y(t), t] = 
d t Y(t), was therefore 



Y(t) = Y(0) exp (iD* 



(9) 



In addition to controlling v y , the variables r and A T were 
kept at zero along the grips, effectively removing any 
shear there. 

At the free boundaries, which comprised X(y, t) in the 
necking simulations and X°(y, t) and X 1 (y, t) in the oth- 
ers, the normal stress a nn and shear stress a„t were set to 
zero, while the tangential stress att was left untouched. 

The simulations were carried out on a uniform rectan- 
gular grid using finite-difference approximations for the 
spatial derivatives. The fields were integrated forward 
in time using the mid-point method, also known as the 
second-order Runga-Kutta scheme, as well as an adaptive 
time-step technique called step-doubling [25| . 

When implementing the STZ equations of motion, a 
numerical viscosity rjV 2 v(() was added to the velocity 
variables on the mapped grid. This was needed be- 
cause the first-order derivatives were discretized using a 
central-difference algorithm; with no higher-order spatial 
derivatives, the grid-points had a tendency to "decouple" 
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into two sets of nodes like the black and white squares 
on a chess board. Note that since the Laplacian was ap- 
plied to the velocities on the mapped grid (that is, on 
the unit square) and not to the usual velocities, the dif- 
fusion term would only be a "proper" viscosity where the 
material was undeformed. 

Since the added viscosity was only meant as a numeri- 
cal tool, it was important to monitor any changes it made 
to the physical results. The influence of the viscosity was 
minimized by performing successive runs with smaller 
and smaller values of r\. If the numerical viscosity was 
chosen too high, the results would look nice and smooth, 
but when compared to simulations with lower values of 77 
it became clear that the added viscosity was smearing out 
sharp features such as shear bands. On the other hand, 
when too small values of 77 were used, large gradients 
caused by the "decoupling" would make the simulations 
more unstable. This was particularly pronounced at the 
border (where values were calculated through extrapola- 
tion) and in larger grids. 

The simulations were most sensitive to the value of 77 
at low strain rates, basically because there was more time 
(per strain) to dampen out the velocities. When running 
as slow as D tot = 10~ 5 on the larger (more unstable) 
grids used in Section there were no values of 77 that 
gave satisfactory results; even the smallest stable value of 
77 would be too large, resulting in unacceptable amounts 
of artificial dissipation and influencing the behavior too 
much. Even for the smaller grids in the necking simula- 
tions of Section llVl the slowest strain rate D tot = 10~ s 
was hard to accommodate. In the end, it was found that 
7/ = 0.02 (which in fact turned out to be a good choice 
for all the simulations) enabled the slowest simulation to 
reach 5% strain while having only a small impact on the 
physical results. The influence due to the exact choice of 
77 for higher strain rates was usually imperceptible. 



IV. NECKING 

Irregular geometry in a two-dimensional, externally 
loaded sample of material can cause stress localizations 
and shear banding. Compared to a zero-dimensional 
model, two spatial dimensions add elastic deformation 
and inertia. As long as the elastic strains are small, 
the speed of sound high, and the rate of total deforma- 
tion low, the difference between a zero-dimensional and 
a uniform two-dimensional system is minimal. Generally 
speaking, metallic glasses have two modes of deforma- 
tion: homogeneous and inhomogeneous 0, 126L WA I2H 
I29I ] . A transition from the former to the latter will occur 
when a material is strained and the stress localizations 
from the irregular geometry start causing shear bands. 
The simulations described in this and the next section il- 
lustrates this transition. Some earlier work can be found 
in 0. 



A. Simulations 

One way of exploring the effect of a non-trivial ge- 
ometry is to look at the dynamics of a material while 
it is necking. A series of simulations were run where 
a rectangular 2x8 piece of material was elongated in 
the y-direction at a constant strain rate.|47| The ma- 
terial was slightly indented in the middle, where the 
width was reduced by 1%. Specifically, the right-hand 
boundary was given by X(y,t = 0) = 1 — 6(y), where 
S(y) = 0.01exp[-ln2(y/0.1) 2 ]. This perturbation to 
the geometry was added to break the symmetry and to 
help trigger any potential instabilities. The material (or 
rather, a quarter of it) was mapped onto an 11 x 41 reg- 
ular grid, using the transformation of variables described 
in Section IIHI The size of the grid was verified to be 
large enough to support the desired accuracy, as a com- 
parison with some test runs on a larger 21x81 grid gave 
almost identical results. The shear modulus was set to 
fi = 100 and the den sity p = 1; since the speed of sound 
is proportional to \fjijp~, a lower value of fi would have 
jeopardized the assumption that the elastic deformations 
were instantaneous. Conversely, an increased value of \x 
would demand notably more resources computationally. 
All the necking simulations had Ao = 1, and the numer- 
ical viscosity was chosen to be 77 = 0.02, based on the 
discussion in Section ITTT1 

Fig. |3 shows the outlines of material samples in two 
separate simulations. The solid outlines were the initial 
configurations; notice that the indentations are hardly 
visible. The dashed outlines represent the configurations 
after the systems had been strained. The simulation on 
the left was strained an order of magnitude faster than 
the one on the right. While the former had a smooth 
neck, the latter had a more irregular boundary due to 
shear bands. 

In general, shear bands would appear more readily 
when the material was strained at a lower rate. The left 
part of Fig.0]shows the plastic dissipation Q as given by 
Eq. © in a simulation with D tot = 10~ 4 , e = 0.03, and 
Ao = 1 after it had strained roughly 6%. The right-hand 
side displays, for the same simulation, the velocity field 
with the uniform part i> U niform = (0,yD tot ) subtracted 
away. These figures show that the material no longer 
was deforming uniformly. Stress concentrations radiating 
out from each notch at roughly 45° increased the plas- 
tic dissipation, and thus the plastic deformation, along 
these shear bands. As can be seen in the velocity plot, 
the whole central piece of the neck became narrower, and 
this is what caused the irregular boundaries seen on the 
right-hand side of Fig. |3| 

In a uniform STZ material, the rate of deformation 
D pl increases drastically when the deviatoric stress |s| be- 
comes larger than unity 0,H3,|^|. In steady state, any 
area with \s\ > 1 can be considered to be flowing plasti- 
cally, while regions with \s\ < 1 are jammed. This applies 
to the tensorial version of the theory, too, when measur- 
ing the deviatoric stress by the invariant s — y/s 2 + t 2 , 



FIG. 3: Examples of the geometry in two numerical experi- 
ments. The solid outlines are the initial configurations — the 
initial indentations are hardly visible. The dashed outlines 
represent the geometry after the left and right samples were 
strained 24% and 10%, respectively. The strain rate in the 
left simulation was D tot = 10 -3 , ten times higher than on 
the right. In both simulations fi — 100, V2 = 0.5, eo = 0.03, 
Aq = 1, and the initial size of the samples were 2x8. 



FIG. 4: The plots in this figure are taken from the same 
simulation as the outline on the right-hand side of Fig. |3] 
but this time at roughly 6% strain. The density plot on the 
left shows the plastic dissipation Q as given by Eq. The 
figure on the right displays the velocity field after the uniform 
part ^uniform = {0,yD tot ) had been subtracted away. (D tot = 
10" 4 , e = 0.03, A = 1) 



also known as the maximum shear stress. Fig. 00 shows a 
series of snapshots from a necking simulation where the 
interior has been shaded according to the values of s; the 
darker the shade, the higher the value (notice that the 
strain is plotted against the total stress a yy at the grip, 
and since a xx = a xy = there, one has that a yy = 2s). 
At low strains the stress was uniform, but as the ma- 
terial necked the stress concentrated in the center. The 
black lines that appear in the last four snapshots mark 
the boundary between s < 1 and s > 1. The steady- 
state solutions suggest that the area between the lines 
were flowing plastically, while the regions on either end, 
outside the lines, were more or less jammed. As the neck 
grew more pronounced, the plastically flowing area would 
shrink, making a smaller and smaller area responsible for 
accommodating the ever increasing global strain. In ad- 
dition, the jammed areas would relax (notice how the 
ends grew lighter in the final snapshots) and contract 
along the strained y-axis while releasing stored energy. 
This decrease in strain at the ends had to be compen- 
sated by increasing the strain even further in the middle. 

The graph below the snapshots plots the true stress 
cfyy at the grips as a function of the total true strain e y ° y . 
The true stress (at the grip) is defined as the force divided 



by the current width of the material, both measured at 
the grip. The true, or logarithmic, total strain is defined 
as £ y °y — ]n[Ly(t)/Ly(0)], where L y (t) is the length of 
the material at time t. The snapshots were marked on 
the graph as circles. The first snapshot was taken while 
the material was in the "elastic phase" , while the second 
was taken just as the material was about to yield. This 
explains why the two first snapshots seem so uniform. 
Note that since the plotted stress was measured at the 
grips, it corresponds to the shading at the ends of the 
material. 

To illustrate how the parameters affect the material 
behavior, Figs. El and [S] show further plots where the 
true stress a yy at the grip was plotted against the true 
total strain e y ° y , comparing curves for various values of 
eo and D tot . In Fig. all the simulations were strained 
at a rate of D tot = 10~ 3 while eo, the amount of strain 
caused by flipping STZs, was varied. In simulations with 
a higher eo, the material needed to be strained further 
before reaching the yield stress. This is because more of 
the work done on the system was stored as plastic energy; 
it was only the elastic deformations that contributed to 
the rising stress. In fact, the slope of the stress-strain 
curves in the "elastic" regime, which corresponds to the 
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FIG. 5: The graph shows the true stress a yy at the grips 
as a function of the total strain e y ° y . Snapshots of the sys- 
tem were taken at six different strains, and the interior was 
shaded according to the value of the maximum shear stress 
s = Vs 2 + t 2 . The black lines in the four last snapshots 
show where s = 1, suggesting that the area between the lines 
(s > 1) was flowing plastically while the areas at the ends 
(s < 1) were jammed. (D tot = 1CT 3 , e = 0.03, A = 1) 



effective shear modulus, has been calculated in Af- 
ter reaching the yield stress, the systems with a higher eo 
would see a faster relaxation of the stress. This is proba- 
bly because a high eo allowed for a higher plastic rate of 
deformation at the neck, where the elastic stored energy 
was released through plastic dissipation. Interestingly 
enough, a high eo meant that most of the stored energy 
was plastic; that is, the energy was stored in flipped STZs 
rather than elastic strain (the amount of stored plastic 
energy cannot easily be deduced from the stress-strain 
curves). 

The effect of varying the total strain rate is shown in 
Fig. [7] The strain rate had little effect on the strain at 
which the yield stress was reached (although it did influ- 
ence the initial stress in these simulations; see below). 




\Y. i I i I I I I I I I 

0.1 0.2 0.3 0.4 0.5 

m . . tot . 

True strain (e ) 

yy 



FIG. 6: The true stress a yy at the grip plotted against the 
true total strain e y ° y for selected values of eo. As eo was in- 
creased, it took longer (in terms of strain) for the material 
to reach the yield stress, since less of the deformation was 
stored as elastic energy (and more as plastic). Once the yield 
stress was reached, a higher value of eo allowed for more plas- 
tic deformation and thus a quicker relaxation of the stress. 
(D tot = I0" 3 , fj, = 100, vi = 0.5, A = 1) 



After the system reached yield stress, the simulations 
with high strain rates took a lot longer (in terms of strain) 
before they would neck. The higher strain rate meant a 
steady-state stress even further above unity (the "yield 
stress" ) , which in turn meant that small geometric inho- 
mogeneities would not be able to separate the stress into 
regions of s > 1 and s < 1. For the low strain rate simu- 
lations, the stress localization was more pronounced, and 
for the lowest rate, the shear bands were so sharp that 
the numerics was not able to strain the material beyond 
5%. 

Fig. [5] compares stress-strain curves of simulations 
where both eo and D tot were varied. These curves show 
clearly that the strain rate primarily changed the time 
between yield and necking, while eo controlled the dy- 
namics of both the stress increase before yield and the 
relaxation after. 

On a side note, all the simulations were started at the 
given (or final) strain rate rather than being ramped up 
from zero. This was especially important for the high 
strain rate simulations. If the material had been started 
from rest, the details of how the strain rate was ramped 
up would greatly affect the outcome of the simulation. If 
ramped up too slowly, the material would yield before the 
final strain rate could be achieved. If the strain rate was 
turned up too quickly, the speed of sound would no longer 
be negligible and inhomogencitics would form at the grips 
due to large stress buildups. To minimize transient effects 
from the non-zero initial strain rate (such as standing 
elastic waves in the material) , initial values for the other 
variables in the simulation were calculated in order to 
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FIG. 7: The true stress a yy at the grip plotted against the true FIQ g . Thg tme stresg ^ at tfafi grip plotted againgt the 



total strain e yy for selected values of D . The value of D 
had little effect on the "rate" , in terms of strain, at which the 
stress reached yield; that is to say, all the simulations in this 
graph yielded approximately at the same strain (the reason 
the D tot = 10 -2 curve starts at a higher stress is explained in 
the text). After yield, though, the high strain rate simulations 
were strained much further before they relaxed (see framed 
plot). (e = 0.03, fi = 100, v 2 = 0.5, A = 1) 



true total strain e % yy . Again, a higher eo caused a slower rise 
to yield but a faster relaxation after. The strain rate did not 
affect the rise to yield, but a faster rate resulted in longer 
plateaus of "uniform steady-state plastic deformation" before 
the indented geometry caused a necking instability, (n — 100, 
u 2 = 0.5, A = 1) 



start the system as close as possible to a steady state. 
This is why the curve with a strain rate of D tot = 10~ 2 
in Fig.[7|starts with a somewhat higher stress. This is not 
a numerical effect; it would be true for real experiments 
as well. See Fig. [5] and the accompanying text for some 
further information on the effects of a non-zero initial 
strain rate. 



B. Discussion 

In the previous section, a series of necking simulations 
with different values of the parameters eo and D tot was 
described. Trends seen in these simulations are discussed 
in more depth in the text below. 

Fig-Elis a "phase diagram" of the dynamics of the neck- 
ing simulations, mapping out the behavior as a function 
of D tot , the total strain rate, and eo, the strain due to 
flipped STZs (Spaepen uses a similar diagram, which he 
calls a "deformation map" , when he outlines the behav- 
ior of amorphous metals in a graph of the temperature 
versus the stress |ll],|2£g). It tries to capture the essence 
of the behavior described in the simulations above. 

Fig- EH shows how the simulations discussed in the pre- 
vious section fit into the "phase diagram" of Fig. [5] Since 
there was more room in this graph, an additional curve 
and shading was added to show what parameter values 
would cause the initial stress to lie between so = 0.1 and 
so = 1, the latter being the yield stress. 

Returning to Fig. the following is an attempt to 
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FIG. 9: A "phase diagram" for the dynamics of constant 
strain rate necking simulations in the (D tot , eo) parameter 
space when A = 1, fi = 100, and V2 = 0.5. Simulations with a 
lower total strain rate had more pronounced shear bands. An 
increased eo resulted in more plastic flow, which allowed for 
sharper necks (in both time and space) . The three curves run- 
ning down the center indicate various values (at steady state) 
of the ratio of stored plastic to stored elastic energy, ip pl /%/j el . 
If parameters were chosen from the shaded area at the top, the 
initial stress would be higher than the yield stress; thus, the 
strain rate would have to be ramped up from a lower value 
(perhaps from rest), and the details of how the strain rate 
was ramped up would affect the outcome of the simulation. 
The simulations in this paper deliberately avoided parameters 
from this region. 
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FIG. 10: This is the same "phase diagram" as in Fig. [5] with 
simulations marked as open and filled circles. The latter are 
the simulations that were presented in Section II V Al The 
shaded area added since Fig. [^represents the parameter space 
where the initial stress would lie between s = 0.1 and so = 1. 
(/j. = 100, u 2 = 0.5, A = 1) 



explain the trends seen there. First, the simulations 
with higher strain rates had less pronounced shear bands. 
Since energy was supplied at a faster rate, the material 
also had to dissipate this energy if it were to stay in 
steady state. In the slow strain rate simulations, it seems 
that the rate was slow enough that all the dissipated en- 
ergy could exit through the shear bands; in other words, 
only the shear bands had stresses above the yield stress. 
In the simulations with higher strain rates the input of 
energy per unit time was higher, resulting in energy dis- 
sipation throughout the material. This forced the stress 
to exceed the yield stress everywhere. Remember that if 
s > 1 then A — > 1/s, while if 5 < 1 then A — > s. Since 
the plastic flow was proportional to eo(s — A) (Ao = 1 
in these simulations), the slowly strained materials saw a 
large difference between the amount of plastic flow in the 
shear bands and the rest of the material. The stresses 
would be above the yield stress in the shear bands and 
below it elsewhere. Notice how the inhomogeneous ge- 
ometry played a vital role by breaking the symmetry, 
concentrating the stresses in certain areas, and thus al- 
lowing the shear bands to form. 

In contrast, the value of the strain rate had little effect 
on the simulations before the stress reached the yield 
stress. As long as A was low, T was negligible, and few 
STZs were created or annihilated. This implies that the 
deformation before yield was almost reversible, and that 
the plastic strain was proportional to A, the fraction of 
flipped STZs. Thus the material would always reach the 
yield stress at the same total (that is, elastic and plastic) 
strain. 

Second, when the parameter eo was given high val- 
ues, the material would have to be strained further be- 



fore reaching yield stress. This parameter controlled the 
amount of strain due to flipping STZs, or alternatively, 
how much of the energy was absorbed in plastic defor- 
mations. By combining Eqs. and Q, one can see 

that a uniform system in steady state with s w 1 has 
^pi/^ei _ eQ ^ + y ^ gjj. when this 

expression was less 
than unity, most of the strain energy was stored as elas- 
tic strain. This would correspond to eo < 0.007 in the 
above simulations, since /i = 100 and V2 = 0.5. As eo was 
increased, the STZs were responsible for more of the de- 
formation, and since the stress is only proportional to the 
elastic strain, the stress would grow more slowly for large 
eo. That explains why the materials with the smaller eo 
yielded earlier. 

Observing the post-yield dynamics, a large value of 
eo ultimately resulted in sharper necks. After the stress 
had reached yield, the energy supplied by the work at 
the grips would no longer be stored in the elastic or plas- 
tic deformations (which corresponds to stress and flipped 
STZs, respectively). Rather, the energy was dissipated 
through flipping newly created STZs while annihilating 
the same amount of already- flipped ones. In a homo- 
geneous system the material would just remain in this 
steady state. In the necking simulations, though, the 
localized stresses emanating from the indentations per- 
turbed the system enough to change the flow from homo- 
geneous to inhomogeneous. Especially for higher values 
of eo, the neck would narrow rapidly, increasing the stress 
and plastic flow there until the material would snap in 
two (in practice, the restricted geometry of the numer- 
ical grid would not actually allow the material to split, 
but the neck would become extremely narrow and the 
time-step would decrease to a value which was, for all 
practical purposes, zero). This uncontrolled necking was 
mainly driven by stored energy being released through 
plastic dissipation at the neck. Interestingly enough, it 
was the stored plastic energy that supplied almost all of 
the energy; the stored elastic energy was almost negli- 
gible when eo 3> 1/^(1 + v-i). An interesting question 
is: How much does the quasi-linear approximation af- 
fect the dynamics? In a fully non-linear version of the 
theory, the STZ transition rates would grow small for 
low stresses, thus preventing the STZs from flipping back 
Potentially, this could prevent some of the release 
of the stored plastic energy from the areas further from 
the neck. 

A large value of eo meant that the material needed to 
be strained further before the stress reached the yield 
stress since the plastic flow relaxed some of the stresses. 
Although it took longer for the material to reach the yield 
stress, once it did the neck developed more rapidly. This 
was because the increased plastic flow in the neck allowed 
the ends to relax and release its stored plastic and elastic 
energy faster. The relaxation of the stress at the ends 
added even more strain to the center, contributing to the 
plastic flow at the neck. Finally, as the neck became thin- 
ner, the stress would rise there and increase the already 
high plastic flow even further. Thus the material would 
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neck faster when eo was high, since that would reinforce 
this feedback loop. 

When eo<C 1/^(1 + ^2)7 most of the initial work done 
on the system was stored as elastic energy. After reach- 
ing yield stress, the steady-state stress became large (and 
A = 1/s small) in order to create enough plastic dissipa- 
tion D pl to counter the steady input of energy supplied 
by the work done at the ends of the material. There were 
some faint shear bands where the stresses showed slightly 
elevated values, but the plastic dissipation was almost 
uniform throughout the material when eg was small. 



V. PRE-ANNEALED MATERIALS 

The behavior of a pre-annealcd metallic glass is more 
brittle than both a virgin as-quenched sample and a ma- 
terial that has experienced plastic work. Experiments 
have shown that metallic glasses that have been annealed 
below the glass transition temperature seem more brit- 
tle, show more pronounced strain softenin g, a nd have 
decreased plastic flow jH |H M, El El 013 • Micro- 
scopically, observations showed structural relaxation in 
the form of more closely packed molecules, correspond- 
ing to a macroscopic increase in density. Upon plastic 
deformation the materials were seen to return to their 
as-quenched states, including a decrease in the density 
and a lower packing fraction. Annealing for a longer 
time or with a higher temperature gave more pronounced 
changes in the mechanical properties. It is also worth 
noting that aging the materials over longer periods of 
time at lower temperatures also made the metallic glasses 
more brittle |27j . 

The STZ theory captures the change in behavior due to 
annealing through A, the relative density of STZs. When 
a virgin material is initially quenched from a molten 
state, the atoms have little time to organize into a closely 
packed configuration, leaving a "fluffy" structure and a 
high density of STZs. Annealing the material packs the 
atoms into a tighter configuration, making Ao < 1. Like- 
wise, if an annealed material is subjected to plastic work, 
the value of A would rise as the tightly packed atoms are 
"re-fluffed" , leaving more room for local rearrangements. 

The version of the STZ theory that is explored numeri- 
cally in this paper has no mechanism for thermal creep or 
relaxation, which is reflected in the fact that A > 0. Thus 
the model can only simulate pre-annealed materials, by 
setting the initial value of A low j^. Consequently, a 
material with Ao < 1 and Ao — 1 will from now on be 
referred to as "pre-annealed" and "worked" , respectively, 
although the latter could also represent a material in its 
virgin as-quenched state. 

Initially during a constant strain-rate experiment, a 
small A would suppress plastic deformation since there 
would be very few STZs to flip. As the stress s — ► 1, 
most of the existing STZs would be flipped, and new 
ones would have to be created causing A to grow. The 
lower eo was set, the higher the peak stress and steady- 



state stress would be, and the faster A would grow toward 
one. 

The previous section assumed that A = 1 through- 
out. The current section describes simulations where 
pre-annealed materials were used during loading, imple- 
mented by setting Ao = 0.01. A spatially extended two- 
dimensional simulation allowed A to grow locally, thus 
making it possible for the material to boost the plastic 
flow in areas where the stresses were high. As will be 
seen, this contribution to the inhomogeneous deforma- 
tion was particularly pronounced for high values of eo- 

As in Section IIVI the simulations had a shear mod- 
ulus of /i = 100, a density of p = 1, a Poisson ratio 
of V2 = 0.5, and a numerical viscosity of r\ = 0.02. 
The size of the material was 4x4, although this time 
the whole material was simulated (that is, there were 
no symmetric boundaries). The material was mapped 
onto a grid measuring 33 x 33 nodes. The left and right 
boundaries X°(y, t) and A 1 (y, t) (initially parallel to the 
y-axis) were both allowed to deform. The lower bound- 
ary was held fixed at y = 0, while the upper boundary 
Y(t) was moved at a constant strain rate (see Section ITTT1 
for further details). The initial density of STZs was set 
below unity, to Ao = 0.01. In order to encourage inho- 
mogeneous flow, the free boundaries were made jagged 
by randomly perturbing the width with values up to 1% 
of the total width (in other words, after the perturbation 
X 1 (y,0) - X°(y,0) £ [3.96,4.04]). By setting the seed 
for the random number generator to the same value for 
all the simulations, the initial geometry would always be 
the same. 

The simulations in this section illustrate the dynamics 
of pre-annealed materials in three different ways. First, 
the distribution and average density of STZs are com- 
pared to the dissipated energy and the work done on the 
systems for different values of e . Second, the dynam- 
ics of pre-annealed and worked materials are compared. 
And third, the following question is addressed: Would 
the simulations with higher strain rates experience more 
pronounced shear bands if they were stopped and held at 
a fixed strain for the same amount of time it would take 
a slower simulation to reach that same strain? 



A. Variation in eo 

Fig, im compares the average density of STZs 

A avg =jjAdA (10) 

(where A is the area of the sample) to the plastic dissi- 
pation of energy 

Qsum= I QdA (11) 
J A 

for three different values of eo- All three simulations were 
strained at a rate of D tot = 10~ 4 . Each of the three 
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graphs on the left displays three curves, of which the 
solid represents the rate of work done on the system at 
the grip, 

External = F V y = (X 1 - X°)a y yV y . (12) 

F is the force applied to the end of the material, and a yy , 
v y , X°, and X 1 were all evaluated at y = Y(t). Since the 
upper boundary was kept flat and moved perpendicular 
to its surface, only the y-component of the stress was 
needed, internal is thus the rate at which energy flowed 
into the system. 

The middle curve in all three graphs shows Q S um, 
the rate of energy dissipated due to plastic deformation 
summed over the whole material. The axis on the left 
side displays the values for both this curve and that of 
^external, although they were multiplied by 10 3 to reduce 
clutter (the interval is really [0,0.004]). The axis on the 
right sets the scale for the third curve, A avg , which repre- 
sents the density of STZs averaged over the whole sample. 

The two density plots to the right of each of the graphs 
show the final distributions of A and Q when the systems 
had reached 7% strain. In the density plots for A, the 
shading was scaled so that the interval [0, 1] went from 
white to black. With Q, the gray scale was mapped onto 
the interval [0, 0.003], with white again representing zero. 
In the latter density plots there were occasional points 
that exceeded 0.003 (these were just painted black as 
well), but choosing a larger interval would have erased 
most of the structure seen in the pictures. 

In the simulation where eo = 0.003, A avg quickly rose 
to one while Q sum grew equal to Pextemai once the stress 
reached yield. The latter meant that all the energy pro- 
vided at the grips was dissipated through plastic defor- 
mation, and none was stored, after the material reached 
about 2% strain. In the eo = 0.3 simulations, A avg never 
even reached a value of 0.15, and the plastic rate of dis- 
sipation grew a lot slower. In fact, looking at the density 
plot shows that the density of STZs only saturated along 
one band, and the plot for Q shows that this was where 
most of the plastic dissipation took place as well. The 
fact that A so quickly reached its equilibrium value every- 
where in the material for small eo meant that the initial 
value of Ao = 0.01 had little effect in this case. 

Although the density plots of Q and A both were snap- 
shots, the latter in some ways gave a cumulative look at 
what had happened since the start, since A could not de- 
cay (this is not true in other versions of the STZ theory 
that incorporate thermal relaxation H2,[2i|). In contrast, 
the density plot of Q is more appropriate when examin- 
ing the instantaneous dynamics, since it highlights the 
current rate of plastic deformation in the different areas 
of the material. 

In the simulation with eo = 0.3, the rate at which en- 
ergy flowed into the system, Pexternab started decaying 
towards the end. This is because the stress at the grip 
began to drop due to the increased plastic flow along the 
shear band. The high value of eo allowed the material to 
deform enough in the small area of this band to account 



for all the strain imposed by the grip, thus inducing inho- 
mogeneous flow and perhaps cause the material to break 
at a later time. 



B. Comparing Pre- Annealed and Worked 
Materials 

Fig. 1121 compares a pre-annealed and a worked mate- 
rial, having initial conditions of Ao = 0.01 and Ao = 1, 
respectively. The two curves in the center show the true 
stress-strain curves of these two materials as they were 
strained beyond yield and toward failure (the materials 
would eventually grow extremely thin somewhere in the 
center). The density plots of Q on either side were arbi- 
trarily chosen some time after the materials had yielded, 
with the stress at the grips being the same in the two 
snapshots. Although the two density plots look different, 
their initial geometries were identical to each other and 
to all the other simulations in this section. 

The only comparison shown here is that of the simu- 
lations with eo = 0.3. Although no less important, the 
pre-annealed materials with a lower eo would have A rise 
to one everywhere so quickly that the results were essen- 
tially the same as for the worked materials; this effect 
was seen in Fig. ^2 

A high value of eo allowed the stresses to grow more 
slowly, permitting some areas to reach s > 1 while others 
remained at s < 1 as the material approached yield. A 
would only grow in the resulting shear bands, thus mak- 
ing its average value grow slowly. A slowly growing stress 
(with respect to time) would also occur for small strain 
rates, resulting in a similar shear-banding effect. 

Comparing the stress-strain curves in Fig. 1121 the pre- 
annealed material plateaued before breaking. It almost 
seemed like it was riding along an unstable equilibrium 
before it increased the density of STZs in one location 
and then deformed and broke there. In comparison, the 
worked material would deform substantially more before 
reaching yield stress, and then "ooze" apart (rather than 
"break"). 

It is tempting, although perhaps somewhat specula- 
tive, to compare the necking on the left side of Fig. IT51 
to fracture. Some of the simulations, including this one, 
had extremely concentrated stresses and narrow necks, 
and the deformations behaved in many ways similar to 
brittle cracks. The simulations have shown that there is 
an interplay between the STZs, the geometry, and stress 
concentrations. Naturally, the geometry of the numer- 
ical grid in these simulations was quite restricted, the 
single-valued boundaries forcing the "crack" to run hor- 
izontally rather than following the 45° shear bands as 
seen in experiments |36ll37| . Also, the resolution was too 
low to make any good quantitative arguments. Never- 
theless, this behavior suggests that the STZ description 
might be capturing some of the dynamics that is present 
in fracture (for more information on fracture in brittle 
amorphous materials, see the review by Fineberg and 
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FIG. 11: Three pre-annealed simulations with Ao = 0.01 were run for eo = 0.003, eo = 0.03, and eo = 0.3, all with a strain 
rate of D tot — 10~ 4 . The graphs on the left show internal, the rate of work done at the grips, Qsum, the rate of dissipated 
plastic energy summed over the whole material, and A avg , the average density of STZs. The left vertical axis shows the range 
for Pextcmai and Q sum (the interval is [0,0.004]), while the right vertical axis displays the scale for A avg . The density plots on 
the right show A and Q at 7% strain. 
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FIG. 12: The true stress-strain curves for two materials, one 
with Ao = 0.01 and one with Ao = 1. The density plots show 
the plastic rate of dissipation at some arbitrary time after 
yield, with the stress at the grips being the same in the left 
and right snapshots. (D tot = 10~ 3 , e = 0.3) 



Marder ,38]). 



C. Strain- Rate Dependent Localization 
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FIG. 13: Density plots of A at 7% strain, both for e = 0.03 
and eo = 0.3. The systems were strained at rates of D ° = 
10 -4 and D tot — 10 -3 , and in the latter case the grip was 
held fixed at 7% for the same amount of time that it took the 
slower strain rate to reach this strain. Despite the extra time, 
no pronounced shear bands developed; in fact, once the grip 
was stopped, the material hardly changed at all. 



When comparing the results of the simulations from 
both the current and previous sections, the shear bands 
at a given strain were more pronounced for smaller strain 
rates. Comparing two identical materials that were 
strained to 7% at rates of D tot = 10~ 4 and D tot = 10~ 3 , 
the former would not only be stretched ten times slower, 
it would also have ten times as long to relax and deform. 
Could it be that it was not the difference in speed, but 
the difference in relaxation time that allowed the shear 
bands to form in the slow case? 

Fig. ^5] shows density plots of A at 7% strain where 
the materials had been strained at the rates mentioned 
above, with both eo = 0.03 and eq = 0.3. The system 
that had been strained at the slower rate, D tot = 10~ 4 , 
reached 7% at time t = 700, while the faster system with 
a rate of D tot = 10 -3 was stopped around t = 70 (the 
snapshot was taken at t = 75 to make sure the material 
had come to a complete stop, reaching its full 7%; in 
order to stop the material at this exact strain, it was 
necessary to start slowing it down already at t = 65). 
After stopping the fast system, it was held fixed until 
t = 700, the same amount of time it took the slow system 
to reach 7% strain. 

It turned out that there was practically no change in 
the quickly-strained material after the grips had come to 
a halt (this was the case both for eo = 0.03 and eo = 0.3, 
although the former ended up with a higher A avg ). This 
means that all the deformation took place instantly, not 
allowing any of the stored energy to escape later. 

This does not preclude a situation where a more irreg- 
ular geometry might induce a stress concentration, allow- 
ing both stored plastic and elastic energy to be released, 
perhaps even driving a necking instability; this kind of 
behavior was seen in the previous section. It does mean, 
though, that the shear bands were created without the 



need of instability mechanisms. In fact, the bands of 
STZs were more pronounced when energy flowed into the 
system at a slower rate. 

The previous section speculated that the STZ theory 
could perhaps contribute to explain fracture dynamics. 
The current section seems to contradict this somewhat, 
at least at first glance. From experience, a material 
strained at a higher rate seems more brittle, implying 
that it should have sharper stress concentrations, and 
more localized shear bands. How could such behavior be 
compatible with the results shown above? 

It was briefly mentioned earlier that the thermal re- 
laxation is bein g in corporated into some versions of the 
STZ theory pll23l]. With that mechanism included, A 
can decrease as STZs are annihilated over time, and this 
might suppress the creation of localized bands of STZs 
as the strain rate drops. The model used in this paper 
did not allow A to decrease, which might be interpreted 
as running the simulations close to zero temperature. 



VI. CONCLUSION 

The goal of this paper has been to explore the STZ 
theory in a spatially extended geometry. In addition to 
showing that a two-dimensional implementation is capa- 
ble of describing shear localization, a wide range of pa- 
rameter values were used to expose the different types of 
behavior inherent in the model. The geometry, energy 
flow, and internal state of the material all contributed to 
the rheology through effects such as jamming and plas- 
tic flow, annealing, strain softening, necking, and shear 
banding. 
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Amorphous metals have two modes of deformation: 
Homogeneous and inhomogeneous flow 0, |2(| |3{| . Ho- 
mogeneous flow is typical for low stresses and high tem- 
peratures, and under uniaxial tension a material sample 
will deform uniformly throughout the specimen. In in- 
homogeneous flow, the deformation is usually localized 
in narrow shear bands that run at a 45° angle with re- 
spect to the tensile axis |36|, |37J (some studies show that 
the angle of the shear bands deviate under large isotropic 
pressures 0,0)- The cross-section of the material de- 
creases as sli p, a nd eventually fracture, occurs along these 
shear bands |26l l36l . At low temperatures, the shear- 
localization instability sets in right after yield, making 
the material behave in a seemingly brittle manner; there 
is no hardening due to strain in metallic glasses, although 
physical aging decreases the plastic response [27| . 

Section llVl looked at how the strain rate, the straining- 
capability of STZs (eo), and the geometry affected flow 
and deformation during necking. It was especially strik- 
ing how sharp shear bands and narrow necks were formed 
at low strain rates and high eo, respectively. The former 
of these two trends appeared to contradict experimental 
evidence: An amorphous metal displays increased brittle 
behavior as it is strained at a higher rate [2jJ , while the 
opposite seemed true for the simulated materials. The 
discrepancy probably stems from the lack of tempera- 
ture in the current STZ model. There was no mecha- 
nism once the simulation had started, apart from plastic 
deformation, to annihilate existing STZs. In the exper- 
iments, the non-zero temperature allowed the molecules 
in the most strained areas to relax if given enough time. 
This resulted in more homogeneous flow as the material 
was strained at lower rates. Efforts have been made to 
incorporate thermal relaxation into the model [22, 23] . 

Section Ivl considered pre-annealed amorphous solids, 
starting with a lower initial relative density of STZs Ao. 
Experiments have shown that amorphous metals become 
more brittle when annealed , ev en though no crystalliza- 
tion is detected [H EJ HI IH HU . Rather, the em- 
brittlement correlates with structural relaxation, leaving 
the more closely packed molecules less room to move. In 
the simulations, the lower Ao enhanced the localization 
of the strain, especially for higher values of eo when most 
of the stored energy was in the form of flipped STZs 
(rather than elastic strain). Compared to either virgin 
as-quenched or worked samples, the pre-annealed mate- 
rials behaved in a more brittle manner with narrower 
shear bands and something resembling a cleavage frac- 
ture. The latter observation should be approached with 



caution since the implementation was not designed to 
handle such extreme deformations; the numerical grid 
forced the "crack" to run horizontally through the ma- 
terial, while a diagonal path running at 45° would seem 
more natural |49j. Nevertheless, the similarities between 
the dynamics of the STZ theory at low Ao and empirically 
annealed solids were strong enough to warrant further in- 
vestigation. 

As mentioned, an increased value of eo would result in 
more of the applied work being stored as flipped STZs 
rather than elastic strain. This was particularly appar- 
ent for the virgin materials (Ao = 1), where an increased 
eo meant that the material would reach yield stress at 
a much larger total strain. On the one hand, it is pos- 
sible that the large values of eo, which were needed to 
produce interesting dynamics in the simulations, exag- 
gerated the strain caused by flipped STZs. On the other 
hand, Ao <C 1 would restrict the plastic strain (and thus 
the stored plastic energy) to narrow shear bands, caus- 
ing a mostly elastic behavior before yield and more brittle 
dynamics at failure. It is possible that the sharp bands of 
inhomogeneous flow in experimentally deforming amor- 
phous solids are due to an internal structure correspond- 
ing to low initial values of A; unless the materials were 
quenched extremely rapidly and to very low tempera- 
tures, chances are that some structural relaxation would 
occur. Incorporating the earlier mentioned thermal re- 
laxation into the STZ model could help lower the value 
of A throughout the simulation and yield better agree- 
ment with the experimental results. 
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